Stripe Orientation in an Anisotropic t-J Model 
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The tilt pattern of the CuOa octahedra in the LTT phase of the cuprate superconductors leads 
to planar anisotropies for the exchange coupling and hopping integrals. Here, we show that these 
anisotropics provide a possible structural mechanism for the orientation of stripes. A t^-ty-Jx- 
Jy model thus serves as an effective Hamiltonian to describe stripe formation and orientation in 
LTT-phase cuprates. 
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Early Hartree-Fock calculations found evidence for 
domain- wall formation in doped 2D Hubbard and t-J 
models. In these calculations the domain walls contained 
one hole per unit cell and separated 7r-phase shifted 
antiferromagnetic (AF) regions. Subsequent density- 
matrix-renormalization-group (DMRG) calculations |^] 
also found hole-domain walls separating 7r-phase shifted 
AF regions, but in these calculations the linear filling of 
the horizontal (or vertical) domain walls corresponded to 
one hole per two unit cells of the wall. In these calcula- 
tions, domain- wall formation originates as a compromise 
in the inherent competition between the kinetic and ex- 
change energies which arises when holes are added to 
a Mott antiferromagnetic insulator. In the parameter 
regime where horizontal or vertical stripes formed, these 
four-fold rotationally invarient models did not distinguish 
between the two orientations. Here we wish to discuss a 
possible electronic mechanism for stripe orientation. 

We are motivated by the structural phase transition ||] 
of Lai,Q-xNdoASrxCu04, in which the system goes from 
a low temperature orthorhombic (LTO) to a low temper- 
ature tetragonal (LTT) phase below ~ 70K. Here, as 
illustrated in Fig. la, in the LTT phase the CuOq octa- 
hedra tilt around an axis oriented along the planar Cu-0 
bonds, say along the y-direction. As a consequence, oxy- 
gen atoms on the tilt axis remain in the plane, but in 
the perpendicular x-direction a staggered tilting pattern 
results with oxygen atoms Oa and Of, in Fig. la displaced 
above or below the Cu02 plane, respectively. The x- and 
y-directions are therefore no longer equivalent in contrast 
to the LTO phase, where the tilt axis is rotated by 45 de- 
grees, as shown in Fig. lb. 

The electronic hopping integrals, and thus the antifer- 
romagnetic superexchange, in the Cu02 planes depends 
sensitively on the Cu-O-Cu bond angle 0. In the specific 
buckling pattern of the LTT phase, this bond angle is 
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9y = t: along the tilt axis direction, but is reduced by 
twice the octahedral tilt angle a in the perpendicular di- 
rection, i.e. 9x = TT — 2a. This x-y anisotropy for the 
electronic hopping and superexchange parameters may 
be conveniently translated into an anisotropic t-J model 
Hamiltonian 

- tx J2 + h.c}j 

{i,i-\-x)a 
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+ JyY: (S. + . • S. - ^) . (1) 
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Here, (i, i-\-x) and (i, i-\-y) denote nearest-neighbor sites 
along the x- and y-directions on a square lattice, respec- 
tively, and doubly-occupied sites are explicitly excluded 
from the Hilbert space. 

The magnitude of the anisotropies is easily estimated 
for typical tilt angles of 4°-5° in Lai g-xNdoASrxCuOi 
with X near 1/8 @. When the tilt axis of the LTT phase 
is vertical, as shown in Fig. la, we have 

— = |cos(7r - 2q;)| and — ^ cos'^ {n ~ 2a) (2) 

ty Jy 

It follows that for a tilt angle of order 4°-5°, At/t ^ 1- 
1.5% and A J/ J ~ 2.-3.%. We note that the direction 
with the larger exchange coupling is naturally also the 
direction with the larger hopping amplitude. Choosing 
t = 500meV and with the exchange coupling constant 
J = 15Q0K of undoped La2Cu04 these estimates give 
At =\tx-ty\r^ 60K and AJ = | - J^j ~ AOK. This 
rough estimate for the exchange anisotropy agrees with 
results from quantum chemistry calculations ^ . 

Given this model Hamiltonian, with Jy > Jx and 
ty > tx, one may ask in which direction stripes are ex- 
pected to form. Since Jy > Jx, the exchange energy is 
optimized by orienting the domain walls along the y-axis 



1 



a.) 



LTT 
tilt axis 





LTO 
tilt axis 



FIG. 1. Planar view of the tilt pattern of the CuOq octahe- 
dra in the (a) LTT and (b) LTO phase. In (a) oxygen atoms 
along the vertical bonds remain in the Cu02 plane while in 
the perpendicular direction they move below (Of,) or above 
(Oa) the plane in a staggered pattern, leading to a reduction 
of tx and Jx relative to ty and Jy. 



tude, we analyze the results of a DMRG calculation to 
obtain further insight in this point. 

We have used DMRG techniques to study a 9 x 8 
lattice with periodic boundary conditions in the 8-site 
y-direction and open boundary conditions in the 9-site 
x-direction. Fig. 2a shows a domain which forms when 
4 holes are added for an isotropic Hamiltonian with 
J /t = 0.35. The boundary conditions cause the domain 
to form around the middle of this 8-leg cylinder. Accord- 
ing to the Hellman-Feynman theorem, 
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Therefore, by calculating the change 
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between the expectation value in the 4-hole ground state 
and the undoped ground state for the isotropic case with 
a given value of J/t, we can determine the variation of 
the domain wall energy with respect to small changes in 
Jx near J. The local change for the individual a;-bonds, 
eq. (H), which contribute to d{H)/dJx are shown on the 
horizontal x-bonds in Fig. 2b for J/t = 0.35. Note that 
these contributions decrease as one moves away from the 
domain wall and we find that 
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In a similar manner we find that the variation of the 
domain wall energy with the exchange energy Jy parallel 
to the wall gives 



niJij 



= 1.20/hole . 
Continuing with the kinetic energy terms, 
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so as to minimize the number of broken exchange bonds 
in the direction with the stronger superexchange. Now, 
one might be tempted to argue that since ty > t^, this 
also lowers the kinetic energy of the system. However, 
transverse motion of the domain walls is also known to 
be important ||l|,^,0|, so that an anisotropy in the hop- 
ping with ty > tx can favor a horizonal orientation of the 
stripes. Because AJ and At are comparable in magni- 
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The local kinetic energy of the 4-hole systems associated 
with the domain wall are listed in Fig. 2c. 
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Alternatively, if the LTT tilt axis runs along the x-axis, 
perpendicular to the domain wall, the shift in energy per 
hole is 

AEa^ = -1.20AJ+ 1.43Ai . (10) 

Therefore, if the domain wall is oriented parallel to the 
LTT tilt axis, there is a net energy reduction (relative to 
an orientation perpendicular to the tilt axes) of 

A£: = 0.7A J - 0.02Ai ~ 25X/hole . (11) 

For a section of domain wall containing 4 holes, this 
would be lOOK. The extensive nature of this energy fa- 
vors alignment of the domain wall with the LTT tilt axis. 

Thus we conclude that the AJ anisotropy dominates 
and favors orienting the stripes along the direction of the 
tilt axis of the LTT phase. This is the same orientation as 
suggested from the "structural corrugation" driven ori- 
entation mechanism originally set forth by Tranquada 
et. al. |]. Here, we have simply looked at a particular 
model in which the corrugation manifests itself by giving 
rise to an anisotropic t-J model. 
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FIG. 2. (a) A 9 X 8 lattice with J/t = 0.35 and 4 holes 
which form a site-centered domain wall. This lattice has pe- 
riodic boundary conditions in the y-direction and open ends 
in the x-direction where a weak, 7r-phase shifted staggered 
magnetic field h = zbO.lt, is applied at the open ends. The 
diameter of the circles indicate the hole density and the length 
of the spins, the spin magnitude. The length of the arrows 
on the left and right sides corresponds to (Sz) ~ 0.37. The 
hole density, which is proportional to the diameter of the gray 
circles, is (n^) — 0.155 along the center of the stripe, (b) The 
change in {{Si ■ Sj — jUirij)) between the 4-hole ground state 
with the domain wall and the undoped Heisenberg lattice. In 
the undoped case the weak applied staggered end-field was 
periodic. The value 0.194 corresponds to the axis of the do- 
main wall, and the values are symmetric about this axis, (c) 
The kinetic energy — ^ ((cL^js + h.c.)) in uints of t of the 
4-hoIe system. 

Now, if the tilt axis of the LTT structure runs along 
the y-axis so that it is parallel to the domain wall, then 
Jy = J, Jx — J — A J, ty = t, and = t — At. In this 
case, the shift in energy per hole of the domain wall due 
to the small anisotropy is 

A^ll = -1.9AJ+ 1.41At . (9) 
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